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Abstract: The isotope composition in precipitation has been widely considered as a tracer of monsoon 
activity. Compared with the coastal region, the monsoon margin usually has limited precipitation with 
large fluctuation and is usually sensitive to climate change. The water resource management in the 
monsoon margin should be better planned by understanding the composition of precipitation isotope and 
its influencing factors. In this study, the precipitation samples were collected at five sampling sites (Baiyin 
City, Kongtong District, Maqu County, Wudu District, and Yinchuan City) of the monsoon margin in the 
northwest of China in 2022 to analyze the characteristics of stable hydrogen (8D) and oxygen (6!8O) 
isotopes. We analyzed the impact of meteorological factors (temperature, precipitation, and relative 
humidity) on the composition of precipitation isotope at daily level by regression analysis, utilized the 
Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT)-based backward trajectory model to 
simulate the air mass trajectory of precipitation events, and adopted the potential source contribution 
function (PSCF) and concentration weighted trajectory (CWT) to analyze the water vapor sources. The 
results showed that compared with the global meteoric water line (GMWL), the slope of the local 
meteoric water line (LMWL; 8D=7.348!8O-1.16) was lower, indicating the existence of strong regional 
evaporation in the study area. Temperature significantly contributed to 8'8O value, while relative humidity 
had a significant negative effect on 8'8O value. Through the backward trajectory analysis, we found eight 
primary locations that were responsible for the water vapor sources of precipitation in the study area, of 
which moisture from the Indian Ocean to South China Sea (ITSC) and the western continental (CW) had 
the greatest influence on precipitation in the study area. The hydrogen and oxygen isotopes in 
precipitation are significantly influenced by the sources and transportation paths of air mass. In addition, 
the results of PSCF and CWT analysis showed that the water vapor source areas were primarily distributed 
in the south and northwest direction of the study area. 
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1 Introduction 


As the main input of terrestrial water resources, precipitation is an important part in the 
hydrological cycle and it is of great significance to economic development and ecological 
stability in inland areas, especially in arid areas. Although heavy isotopes of hydrogen (87H, also 
known as 5D) and oxygen (5/80) in precipitation exist with very small absolute amounts, they are 
very sensitive to environmental changes. They can record the information of water transportation, 
condensation, exchange, and other atmospheric process while the air mass is transported from its 
source to the site of precipitation (Craig, 1961; Dansgaard, 1964). As a natural tracer in 
precipitation, the stable 57H and 5!°O isotopes are widely used to study moisture sources (Tan, 
2014; Wang et al., 2017; Fang et al., 2022; Geppert et al., 2022; Sun et al., 2023; Tharammal et 
al., 2023). 

Gimeno et al. (2012) summarized the methods available to identify the moisture sources, 
including the analytical box method, physical tracer (isotope) method, and numerical tracer 
method. Among them, the physical water vapor tracer method, which includes the vertical and 
horizontal atmospheric processes, enables the combination of global climate models and 
Lagrangian models. In recent years, this method has been widely used to analyze the source of 
water vapor by mapping the movement of air mass containing water vapor (Breitenbach et al., 
2010; Chen et al., 2020; Cloux et al., 2021). 

In precipitation, stable &H and 8!8O isotopes are components of water molecules, and their 
composition variations depend on moisture source, meteorological factor, and transportation 
mode (Gou et al., 2022). The regimes controlling isotope fractionations in precipitation have been 
concerned by a large number of researchers (Baker et al., 2015; Steen-Larsen et al., 2015; Tang et 
al., 2017; Espinoza et al., 2021). Due to the large precipitation amount, many studies on the 
influencing mechanism of water vapor source on precipitation isotope have been conducted in 
monsoon area. According to the convective activity, cloud formation, and air trajectory in the 
monsoon area, the researchers examined the influence of various water vapor sources on the 
variation of isotopes in precipitation on different scales (Cai et al., 2018; Lone et al., 2020). For 
example, the decrease in the 5!%O value of precipitation on the Tibetan Plateau, China at the 
beginning of monsoon season (from late May to early June) can be attributed to the changes in 
water sources (Yao et al., 2013). Moreover, convective processes may change the composition of 
water vapor isotopes in the source region, and the transport of water vapor also affects the 
composition of isotopes in precipitation (Lee et al., 2007; Worden et al., 2007; Risi et al., 2008; 
Kurita, 2013; Moore et al., 2014). From a spatial perspective, the temperature effect, amount 
effect, and altitude effect are important factors affecting isotope composition (Salamalikis et al., 
2016; Jiao et al., 2019; Yang et al., 2019; Pant et al., 2021). However, many studies found an 
association between the composition of isotopes in precipitation and the location and/or intensity 
of convection activities in the upper reaches (Zhou et al., 2019). In addition, the falling raindrops 
undergo isotope fractionation when passing through the unsaturated air, namely secondary 
evaporation under the cloud base, which changes the composition of isotopes (Chen et al., 2020). 
The value of deuterium excess (d-excess, which is equal to 5D-85'%O) is related to the 
evaporation rate of precipitation source area, and a high value of d-excess usually indicates low 
relative humidity and intense evaporation (Natali et al., 2022). 

Compared with the coastal region dominated by monsoon, the monsoon margin usually has 
limited precipitation amount with large fluctuation and is usually sensitive to climate change. The 
concept of "monsoon triangle" was first proposed by Li et al. (1988) who pointed out that a 
sensitive triangle area to climate exist in China during the Quaternary period, and the great 
changes in circulation pattern, vegetation, soil, and climate have been found in the "monsoon 
triangle" between glacial and interglacial periods. The western vertex of "monsoon triangle" 
reaches deep into continent, while the other two vertexes are close to sea. Actually, the western 
vertex corresponds to the monsoon margin. After a long-range transport, air mass sources 
reaching the monsoon margin are usually complex, and frequent precipitation extremes can be 
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reported in such an arid climate. In addition, the climate research in the monsoon margin is also 
an important topic, and the findings from modern precipitation isotopes are necessary basis for 
paleoclimate proxy. A systemic and synchronous observation network of precipitation isotopes in 
the monsoon margin not only enhances the understanding of evaporation and condensation in arid 
condition, but also has great significance for ecological protection and high-quality development 
in water scarcity areas. 

The objectives in this study are to: (1) use the Lagrange method to examine the linkage 
between the water vapor sources and precipitation stable isotopes (5!8O and 5D) in the monsoon 
margin in an arid condition; and (2) quantify the influence of different water vapor sources on 
isotopic characteristics. 


2 Data and methodology 
2.1 Study area 


Here five meteorological stations (Fig. 1) were selected from the margin area of East Asian 
monsoon (also known as the margin area of "monsoon triangle" or monsoon margin) as sampling 
sites in this study, including Baiyin City, Kongtong District, Maqu County, and Wudu District in 
Gansu Province, as well as Yinchuan City in Ningxia Hui Autonomous Region, China. The 
highest elevation is in Maqu County (3471 m a.s.l.), and the lowest is 1010 m a.s.l. in Yinchuan 
City (Table 1). The study area has a dry continental monsoon climate, with uneven distribution of 
annual precipitation and intense evaporation (Kong et al., 2023). The average annual temperature 
in the study area is 9.29°C, and the average annual precipitation is 394.66 mm, with more rains in 
summer and autumn, and less rain and snow in spring and winter. This network corresponds to 
the western vertex of "monsoon triangle" of China (Li et al., 1988). 


Fig. 1 Location of sampling sites at the margin area of the East Asian monsoon 


Table 1 Climate characteristics of the five sampling sites from 1991 to 2020 
Sampling site Latitude Longitude Altitude (m) Temperature (°C) Precipitation (mm) Relative humidity (%) 


Baiyin City 36°55'N 104°18'E 1719 9.03 210.40 56.00 
Kongtong District 35°55'N 106°67'E 1364 9.70 486.40 51.00 
Maqu County 34°00'N 102°08'E 3471 2.30 612.30 63.00 
Wudu District 33°00'N 104°92'E 1116 15.30 470.40 61.00 
Yinchuan City 38°48'N 106°22'E 1010 10.10 193.80 52.00 


Note: The latitude and longitude data are the geographical coordinates of the meteorological stations that were used as the sampling 
sites by this study. 
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2.2 Isotope data 


A funneled polyethylene bottle with a capacity of 1 L was placed at the sampling sites as a 
precipitation collecting device. To avoid evaporation, we collected the precipitation samples after 
each precipitation event. In this study, precipitation event was defined as more than 1 mL of water 
was collected in the bottle. After collection, the collected water samples were transferred to sealed 
50 mL polyethylene bottles, and then kept in a cool place. For solid precipitation, it was placed 
into a zip-locked bag at room temperature after collection, and then transferred to a polyethylene 
bottle after melting. The basic information of meteorological elements (including temperature, 
humidity, wind direction, wind speed, etc.) was recorded when precipitation event occurred. A 
total of 236 precipitation samples were collected at the five sampling sites in the whole year of 
2022. The water samples were measured in the Stable Isotope Laboratory, College of Geography 
and Environmental Science, Northwest Normal University, China The DLT-100 liquid water 
isotope analyzer (Los Gatos Research, San Jose, the USA) was used to measure the values of 'H, 
°H, 7H, '°O, 5!70, and!8O. The analysis error of 5D value was within +0.6%o and that of the 5!°O 
value was within +0.2%o in this study. The values of ôD and 5!8O were expressed as 
one-thousandth of the difference relative to the Vienna Standard Mean Ocean Water (VSMOW) 
(Zhang et al., 2023): 


Rsample 


ô 


Rstandard : i () 
where 6 is one-thousandth of the difference relative to VSMOW; Rsample is the ratio of heavy to 
light isotopes in water samples, i.e., D/H and !8O/'°O for SD and 5!80, respectively; and Rstandard is 
the corresponding isotopic ratio in VSMOW. 

Surface meteorological data including temperature, wind direction, wind speed, and 
precipitation were obtained from the China Meteorological Data Network (http://data.cma.cn/dat). 
The upper air meteorological data obtained from the National Centers for Environmental 
Prediction (http://ready.arl.noaa.gov/archives.php) were used for air mass backward trajectory 
analysis as well as potential source area analysis. 


2.3 Methodology 


In this study, Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model was 
used to simulate the air mass trajectory of precipitation events. The potential source contribution 
function (PSCF) and concentration weighted trajectory (CWT) were used to analyze the water 
vapor sources of precipitation air mass. 


2.3.1 HYSPLIT model 


The HYSPLIT model developed by the National Oceanic and Atmospheric Administration 
(NOAA) and the Australian Bureau of Meteorology was used to simulate patterns of airflow 
movement, deposition, and dispersion. It has been widely used to track sources and paths of 
precipitation air mass. In the model, the Lagrangian method is used for advection and diffusion, 
and Euler's method is used for concentration calculations (Stohl et al., 2005). MeteoInfo 
backward trajectory clustering analysis based on the HYSPLIT model was used in the water vapor 
source analysis (Wang, 2014). The Global Data Assimilation System (GDAS) data were imported 
into MeteolInfo software (Chinese Academy of Meteorological Sciences, Beijing, China). We 
conducted backward trajectory simulation for precipitation events in the study area, and then 
clustered the trajectories obtained. 


2.3.2 PSCF 


The PSCF analysis is a method used to identify possible source locations based on clustering 
results of backward trajectory simulation. Considering the environmental significance of d-excess 
(Natali et al., 2022; Chen et al., 2023; Kong et al., 2023; Pérez-Alarcon et al., 2023; Tao et al., 
2023), we applied the average d-excess value as a threshold to analyze the precipitation potential 
source area in this study. That is, the trajectories with d-excess values higher than the average 
d-excess value were counted as potential trajectories; otherwise, they were not counted as 
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potential trajectories (Shi et al., 2010; Xie et al., 2022). The PSCF value is the ratio of number of 
trajectories that pass through a target grid (m) to the number of all trajectories (n), with the 
following equation: 
PSCF=”. (2) 
n 
If the number of trajectories passing through the target grid (m) is too small, the probability 
function results will deviate from the actual results. To ensure the accuracy of PSCF value, we 
introduced a weighting factor w defined in Zeng and Hopke (1989) as follows. 


WPSCF = peek. (3) 
w 


where wpscr is the weight of PSCF value. 
2.3.3 CWT 


Besides PSCF analysis, we also adopted CWT analysis to identify the potential source location of 
trajectory in this study (Hsu et al., 2003). The spatial distributions of trajectory and corresponding 
isotopic characteristics were considered in this method (Zannoni et al., 2019; Zhang et al., 2023). 
The formula for calculating CWT is as follows: 
Dae CkTijk 
Cy = = (4) 


n 
__, Cijk 
i=l 


where Cy is the weighted average of d-excess value in the grid (i, J) (%o/m°); k is the trajectory; n 
is the total number of trajectories; Cx is the d-excess value when the trajectory k travels in the grid 
(i, j); and tj is the duration that the trajectory k stays in the grid (i, j). High CWT value (i.e., Cj) 
indicates a large potential contribution of a specific location. If the numbers of trajectory endpoint 
passing through the grid are too small, it will lead to biased results. To ensure the accuracy of 
CWT value, we used the same weighting factor w as PSCF (Zeng and Hopke, 1989), and the 
formula is as follows: 

wewt = wx CWT. (5) 
where wcwr is the weight of CWT value. 


3 Results 


3.1 Isotope characteristics 


3.1.1 Spatial and temporal variations of hydrogen and oxygen isotopes 


The summer monsoon transports water vapor from the tropical or subtropical oceans to inland 
areas (Ma et al., 2022). The monsoon moisture usually starts in mid-May, gradually advances 
northward, and ends completely in mid-October (Bian et al., 2017). During June—September 
2022, the lowest monthly precipitation was 2.60 mm (September in Baiyin City) and the highest 
was 246.90 mm (August in Kongtong District). July and August were periods with relatively 
frequent precipitation events, and monsoon circulation brought a lot of water vapor to the 
monsoon margin. The temperature fluctuated between 3.50°C and 29.50°C during June- 
September 2022, and the average was 18.20°C. Among the sampling sites, the temperature in 
Maqu County was relatively low, with a fluctuation range of 3.50°C—17.10°C and an average of 
11.70°C. This may be mainly because Maqu County has a high altitude with a plateau climate, 
whose characteristic is long cold winters and short cool summers. The 5D value in June- 
September 2022 ranged from —145.35%o0 to 47.18%, with an average value of —37.79%o, and the 
5/80 value ranged from —19.34%o to 7.42%o, with an average value of —5.31%o (Table 2). These 
results are within the range of the statistical results of stable hydrogen and oxygen isotopes 
(—263.20%o-59.00%o for 5D, and —33.40%o—8.50%o for 5'8O) of precipitation in the arid region of 
Northwest China (Zeng et al., 2020). 
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Table 2 Statistic values of stable hydrogen (6D), oxygen (8!8O), and deuterium excess (d-excess) for each 
sampling site during June—September 2022 


SD (%o) 5'8O (%o) d-excess (%o) 
Sampling site - - : 
Max Min Mean SD Max Min Mean SD Max Min Mean SD 

Baiyin City 34.53 -120.02 -17.65 32.84 4.69 -16.02 -2.22 449 16.01 -14.37 0.11 8.25 
Kongtong District 28.81 -118.99 -40.87 31.95 3.12 -16.37 -640 4.18 21.17 -8.88 6.22 6.92 
Maqu County -3.90 -120.68 -71.06 17.36 -2.21 -15.02 -9.38 3.14 1845 -27.28 3.94 13.85 
Wudu District 47.18 —84.22 -28.13 29.81 7.42 -11.37 -3.79 4.05 20.67 -25.25 2.16 9.81 
Yinchuan City 35.22 -145.35 -31.25 38.75 5.96 -19.34 -4.76 5.23 22.07 -20.11 7.21 9.97 
Total 47.18 -145.35 -37.79 34.94 7.42 -19.34 -5.31 4.79 19.67 -19.18 3.93 11.27 


Note: Max, maximum value; min, minimum value; SD, standard deviation. 


3.1.2 Local meteoric water line (LMWL) 


Craig (1961) found a linear correlation between 5D and 5!%O isotopes in natural water samples 
with the equation of 5D=8.005!8O+10.00, also known as the global meteoric water line (GMWL). 
Rozanski et al. (1993) updated the line using the global precipitation isotope data as 
5D=8.175!80+10.56. According to the 1815 months of observed precipitation isotope data in 
China, the Chinese meteoric water line is 5D=7.805'80+8.70 (Wang et al., 2022), and the 
nationwide grid isotope data with finer spatial resolution (the number of grids is 431,136) show 
an equation as 5D=7.408!8O+5.50 (Wang et al., 2022), both the lines are generally similar to the 
global mean. 

Due to the differences in climatic conditions, the isotopic composition of precipitation spatially 
varies. Strong evaporation causes the water line equation to deviate from the slope of 8 of 
GMWL. The higher the temperature is, the more intense the sub-cloud evaporation is, and the 
smaller the slope and intercept are (Craig, 1961). The LMWL of the summer half year (from 
March to September) was 5D=6.855!8O-2.11, and that of the winter half year (from October to 
February) was 5D=6.665!8O-7.48 (Fig. 2a). Some discrete points in Figure 2 were distributed on 
the right side of meteoric water line, indicating that there was a strong sub-cloud secondary 
evaporation in the study area. The slope of LMWL in summer was higher than that in winter, 
indicating more precipitation and higher relative humidity in summer than in winter in the study 
area. As shown in Figure 2b, the slope of LMWL in Kongtong District was 8.16, followed by 
Baiyin City (8.12), Wudu District (7.43), Yinchuan City (6.79), and Maqu County (4.68); and the 
intercept of LMWL in Kongtong District was 10.52, followed by Yinchuan City (0.75), Wudu 
District (-1.27), Baiyin City (—2.46), and Maqu County (—-27.75). The equation of LMWL in the 
whole study area was 5D=7.345!8O-1.16 (Fig. 2c), which was lower than the global (Rozanski et 
al., 1993) and national (Wang et al., 2022) levels. This is mainly because of the low humidity and 
strong evaporation in arid area, which leads to the enrichment of stable 8D and 8!80O isotopes in 
precipitation. 

3.1.3 Environmental effect 


Meteorological factors such as temperature, precipitation, and relative humidity are considered as 
the main environmental factors affecting the composition of precipitation isotope. Therefore, in 
this study, we discussed the relationship between 5!%O value in precipitation and meteorological 
factors at daily level. The results showed that there was a correlation between 5'8O value and 
temperature (5!8O=0.24Tem—9.59, R?=0.15, P<0.001; Tem, temperature). The temperature effect 
for each sampling site was weaker than in the whole study area. There was a significant positive 
correlation between 5!%O value and temperature in Kongtong District (R°=0.13, P<0.050). Maqu 
County is located in a continental alpine region with a low temperature throughout the year and a 
long cold season; in contrast, Wudu District has a high temperature throughout the year with 
sufficient precipitation and high relative humidity; we didn't find significant temperature effect at 
these two sampling sites. There was a significant negative correlation between 5!°O value and 
relative humidity in the margin area of East Asian monsoon (8!8O= —0.16RH+5.09, R?=0.18, 
P<0.001; RH, relative humidity). A negative correlation was also found between 5'8O value and 
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precipitation (8'8O= —0.05Pre—5.82, R°=0.01, P=0.104; Pre, precipitation), although the 
correlation between 5'8O value and precipitation at most sampling sites was not statistically 
significant at 0.05 level. 
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Fig. 2 Relationship between hydrogen (SD) value and oxygen (8!8O) value of precipitation in the study area in 
2022. (a), relationship between 5D value and 5!80 value of precipitation during summer half year and winter half 
year of 2022; (b), relationship between ôD value and 8!8O value of precipitation at each sampling site in 2022; 
(c), relationship between 5D value and 5!8O value of precipitation in the whole study area in 2022. LMWL, local 
meteoric water line; GMWL, global meteoric water line. 


3.2 Backward trajectory analysis 


The HYSPLIT-based backward trajectory model is considered as an effective method to simulate 
airflow movement. The backward duration for each trajectory was calculated as 10 d (240 h) 
before precipitation events, which is nearly the upper limit of the time that water vapor stays in air 
(Trenberth, 1998; Numaguti, 1999). The trajectory was calculated at 6-h intervals on the day of 
precipitation event occurred, that is, 00:00, 06:00, 12:00, and 18:00 (Beijing Time). We chose 
1500 m above ground level as the start height for each sampling site. The Baiyin City and Wudu 
District were mainly affected by continental water vapor, and the northwest continental water 
vapor had the greatest influence for Baiyin City, accounting for 49.93% of the total water vapor 
(Fig. 3). The air mass in Kongtong District mainly came from the northwestern (37.36%), 
southeastern (23.70%), and western (23.28%) continent, and 15.66% of water vapor came from 
the Black Sea. The Maqu County was mainly affected by the westerlies, including 31.50% from 
the western continent, 28.68% from the Mediterranean Sea, and 13.82% from the Atlantic Ocean. 
In addition, Maqu County was also affected by water vapor from the southeast continent and the 
Indian Ocean. For Wudu District, the largest influence was from the western continent, 
accounting for 48.63% of the total water vapor. The Yinchuan City was mainly affected by water 
vapor from the northwest continent (78.70%), of which 24.32% of water vapor came from the 
direction of Russia. In addition, Yinchuan City was affected by a small part of local water vapor 
cycle (21.30%). 
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Fig. 3 Air mass trajectory (a, c, e, g, and i) and its clustering result (b, d, f, h, and j) at each sampling site in 
2022. Different color lines in the right panel indicate that air mass comes from different directions; of which 
yellow line represents air mass comes from the northwest continent, grey line represents air mass comes from the 
northwest continent in the direction of Russia, red line represents air mass comes from the western continent, blue 
line represents air mass comes from the Black Sea, green line represents air mass comes from the Mediterranean, 
pink line represents air mass comes from the Atlantic Ocean, purple line represents air mass comes from the 
southwest continent, orange line represents air mass comes from the southeast continent, and white line represents 
air mass comes from the proximity continent. Note that the base map used for computerizing this map is from the 
Meteoinfo website (http://www.meteothink.org/index.html#). 
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3.3 Potential source area analysis 


3.3.1 Identification of moisture sources 


To further explore the sources of water vapor in the monsoon margin, we simulated the backward 
trajectories of 236 precipitation events from January to December 2022 in this study. Water 
samples were collected when every precipitation event occurred, and a precipitation event was 
defined as more than 1 mL of water was collected in this study. As the common operation in other 
studies, we took 10 d as the duration of trajectory before the occurrence of precipitation event. 
Based on the results of backward trajectory of HYSPLIT model, we divided the moisture sources 
of precipitation events into eight types (Table 3; Fig. 4): continental moisture from northwest 
(CNW), continental moisture from west (CW), continental moisture from proximity region (CP), 
continental moisture from southwest (CSW), continental moisture from southeast (CSE), moisture 
from the Indian Ocean to South China Sea (ITSC), moisture from the Atlantic Ocean (OA), and 
moisture from polar region (PR). 


Table 3 Climatic, 8D, and 6!8O characteristics of each kind of water vapor source for the monsoon margin in 
2022 


Water vapor j Length v (m/s) Tem RH Pre 8D (%o) 5O (%o) d-excess (%o) 
peur’, (km) CC) (%) (mm) Mean SD Mean SD Mean SD 
CNW 33 7133 2.08 15.00 63.58 185.00 -23.24 31.65 -3.39 457 3.85 11.01 

CW 76 6011 2.10 13.20 62.37 330.80 -39.83 37.70 -5.53 5.05 440 12.20 
CP 13 3140 2.02 10.40 79.60 45.00 -46.30 28.14 -6.11 3.66 2.57 12.70 
CSW 10 3152 2.19 9.70 76.50 31.80 -62.38 33.84 -8.22 5.62 3.38 13.86 
CSE 9 3611 1.58 13.90 84.90 97.60  -36.88 22.38 -6.17 2.62 12.49 5.02 
ITSC 79 5078 10.72 15.90 78.20 772.30 -58.28 24.81 -7.96 3.58 5.41 10.64 
OA 10 16,815 2.30 8.60 70.79 37.90 -56.98 42.70 -7.54 5.62 3.32 12.14 
PR 6 13,121 1.88 14.40 63.10 44.70 -9.55 23.41 -144 3.40 1.97 11.69 


Note: n represents the number of precipitation events related to each kind of water vapor source in 2022. v, wind speed; Tem, 
temperature; RH; relative humidity; Pre, precipitation; SD, standard deviation; CNW, continental moisture from northwest; CW, 
continental moisture from west; CP, continental moisture from proximity; CSW, continental moisture from southwest; CSE, continental 
moisture from southeast; ITSC, moisture from the Indian Ocean to South China Sea; OA, moisture from Atlantic Ocean; PR, moisture 
from polar region. 


For CNW, the average length of air mass trajectory was 7133 km. In 2022, 33 precipitation 
events related to CNW occurred, contributing 13.98% of the total 236 precipitation events in the 
monsoon margin, and such precipitation events were mainly concentrated in the summer half year 
(May-September). This part of air mass contributed about 12.00% of the total precipitation. This 
type of air mass originated from the northwest continent such as northern Siberia, Mongolia 
Plateau, and other regions, and extended to the study area; and it affected precipitation in the 
study area with partial local recycled water vapor. The air mass from this direction had low 
relative humidity, and high values of 5D and 8!%O, and brought little precipitation to the study 
area. 

For CW, the average length of air mass trajectory was 6011 km. In 2022, there were 76 
precipitation events related to CW, contributing 32.20% of the total number of precipitation 
events and 21.40% of the total precipitation in the monsoon margin. This type of air mass 
originated from Western Eurasia and extended to the study area by horizontal transport. It can be 
seen that the study area was mainly influenced by westerly air mass, from the results of backward 
trajectory clustering analysis done by this study. Especially in winter, CW contributed a high 
percentage of precipitation and a high relatively d-excess value. 

For CP, the average length of air mass trajectory was 3140 km. In 2022, 13 precipitation events 
were related to CP, contributing 5.51% of the total number of precipitation events and 2.90% of 
the total precipitation in the monsoon margin. This kind of air mass was mainly from the 
continental area near the study area, and the precipitation events related to CP were mainly 
concentrated in the period around the monsoon (May—October). This part of the air mass had high 
relative humidity but produced little precipitation and the lowest d-excess value. 
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Fig. 4 Water vapor sources of precipitation in the monsoon margin in 2022. (a), continental moisture from 
northwest (CNW); (b), continental moisture from west (CW); (c), continental moisture from proximity region 
(CP); (d), continental moisture from southwest (CSW); (e), continental moisture from southeast (CSE); (f), 
moisture from the Indian Ocean to South China Sea (ITSC); (g), moisture from the Atlantic Ocean (OA); (h), 
moisture from polar region (PR). The denser the line, the more water vapor from this direction. Note that the base 
map used for computerizing this map is from the Meteoinfo website (http://www.meteothink.org/index.html#). 


For CSW, the average length of air mass trajectory was 3152 km. In 2022, there were 10 
precipitation events related to CSW, contributing 4.24% of the total 236 precipitation events and 
about 2.10% of the total precipitation in the monsoon margin. This kind of air mass was mainly 
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from the southwest continent near the study area like India and the Bay of Bengal. The 
precipitation was mainly concentrated in October after the summer monsoon. This part of air 
mass had high relative humidity and produced the lowest values of 5D and 8/80. 

For CSE, the average length of air mass trajectory was 3611 km. In 2022, 9 precipitation events 
related to CSE occurred, contributing 3.81% of the total number of precipitation events and about 
6.30% of the total precipitation in the study area. This kind of air mass was mainly from the 
southeast continent of the study area and the Yellow Sea. Precipitation mainly came from the 
southeast monsoon, which happens in June—September each year. This part of air mass had high 
relative humidity, high temperature, and strong evaporation effect, and it produced the highest 
d-excess value. 

For ITSC, the average length of air mass trajectory was 5078 km. In 2022, there were 79 
precipitation events related to ITSC, contributing 33.47% of the total 236 precipitation events and 
about 50.00% of the total precipitation in the study area. This kind of air mass was mainly from 
the Indian Ocean, Arabian Sea, and South China Sea. Precipitation from this kind of air mass 
mainly dropped during June—September. At the beginning of monsoon period, carrying a lot of 
water vapor, air mass from ITSC passed through the southwest of China at a high speed and 
brought abundant precipitation to the study area. It usually had low values of 5D and 5'80. 

For OA, the average length of air mass trajectory was 16,815 km. In 2022, there were 10 
precipitation events related to OA, contributing 4.24% of the total number of precipitation events 
and about 2.50% of the total precipitation in the study area. This type of air mass originated from 
the Atlantic Ocean and Mediterranean Sea and extended to the study area by zonal transport. 
Precipitation related with this kind of air mass mainly dropped during May-July. This kind of air 
mass had a relatively low humidity. However, the long-distance water vapor transportation caused 
the heavy isotopes in water vapor to be gradually consumed, resulting in a lower d-excess value. 

For PR, the average length of air mass trajectory was 13,121 km. In 2022, 6 precipitation 
events related to PR occurred, contributing 2.54% of the total number of precipitation events and 
about 2.90% of the total precipitation in the study area. This type of air mass mainly originated 
from the Arctic Ocean and was transported to the study area through the northwest continent. 
Although it came from ocean, this type of air mass had low relative humidity, due to the long 
travel path. After passing over the Siberia and Mongolian Plateau, this kind of air mass had the 
lowest d-excess value. 


3.3.2 Identification of water vapor source areas by PSCF and CWT 


The PSCF analysis determined the magnitude of the impact of potential water vapor source areas, 
that is, the proportion of trajectories exceeding the threshold value of d-excess for each grid. Then 
CWT analysis was used to identify the magnitude of the weighted concentration values of 
potential evaporation source areas. The water vapor in Baiyin City mainly came from the west 
and southwest of the study area (Fig. 5). The water vapor in Kongtong District mainly came from 
the surrounding area and southwest direction. The precipitation sources in Maqu County were 
mainly distributed in the southwest of the study area. The water vapor in Wudu District mainly 
came from the south part. The water vapor sources of Yinchuan City were mainly distributed in 
the south and southeast. The results of CWT analysis showed that Baiyin City and Wudu District 
were mainly affected by CW and ITSC; while Kongtong District, Maqu County, and Yinchuan 
City were mainly affected by CP, ITSC, and CW. 

3.3.3 Relationship between d-excess value and wind 


The d-excess value responds to the climatic environment of water vapor source (Dansgaard, 
1964). When the moisture source is relatively dry with low humidity and rapid evaporation, the 
d-excess value will be relatively high; conversely, the d-excess value will be relatively low (Malik 
et al., 2022). Figure 6 shows the relationship between d-excess value (interpolated using all 
precipitation events) and wind characteristics (wind speed and wind direction) for all precipitation 
events at the five sampling sites in 2022. High d-excess values in Baiyin City mainly appeared in 
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Fig. 5 Results of potential source contribution function (PSCF; a, c, e, g, and i) and concentration weighted 
trajectory (CWT; b, d, f, h, and j) for the potential water vapor source area analysis of the selected five sampling 
sites in the monsoon margin in 2022. Note that the base map used for computerizing this map is from the 
Meteoinfo website (http://www.meteothink.org/index.html#). 
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Fig. 6 Rose diagrams showing the relationship of deuterium excess (d-excess) value with wind speed (v) and 
direction at each sampling site. (a), Baiyin City; (b), Kongtong District; (c), Maqu County; (d), Wudu District; (e), 
Yinchuan City. 


the precipitation of air mass blown from the south and southwest of China, and a small part of 
precipitation from the southeast of China when the wind speed is lower (<1.50 m/s). High 
d-excess values of Kongtong District mainly appeared in the precipitation of air mass blown from 
the southeast and southwest of China, with a small part from the northwest of China, when the 
wind speed is high (3.00—4.00 m/s). Precipitation was abundant in Maqu County, and high 
d-excess values mainly appeared in the precipitation of air mass originated from the southeast and 
southwest of China, and a small amount of precipitation from the northeast of China. High 
d-excess values of Wudu District mainly appeared in the precipitation of air mass originated from 
the southeast, south, and southwest of China with lower wind speed (<1.50 m/s). High d-excess 
values of Yinchuan City mainly appeared in the precipitation of air mass originated from the 
southeast and southwest of China. 


4 Discussion 


4.1 LWMN in the monsoon margin 


The composition of 8H and &!ŝO isotopes in precipitation air mass varies greatly, and the 
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difference of 57H and 5'8O isotopic composition can reflect the source of water vapor and the 
degree of water vapor evaporation (Claus et al., 2021; Bhat et al., 2022; Lekshmy et al., 2022). In 
this study, the equation of atmospheric water line in the western fringe of "monsoon triangle" was 
derived by the least square method. Compared with the GMWL, the slope of LMWL in the study 
area was relatively low. The results showed that the precipitation process was accompanied by 
secondary evaporation in the study area. In addition, the secondary evaporation was more intense 
in winter half year than in summer half year in the study area, according to the difference of 
d-excess value between these two periods, to be specific, the d-excess value of winter half year 
was lower than that of summer half year due to the lower relative humidity in winter. 


4.2 Identification of water vapor sources 


Precipitation in an area mainly consists of local evaporative water vapor and long-distance 
transported water vapor. Water vapor transported from outside the precipitation area is essential 
for the generation of intense precipitation (Zhao et al., 2023; Xie et al., 2023; Vishwakarma et al., 
2024), especially for the northwest inland areas where air humidity is low, while evaporation is 
intense. Understanding the specific influence of different water vapor sources could better grasp 
the water cycle process in arid area. Zeng et al. (2020) found that westerly water vapor is the main 
water vapor in the northwest, which is consistent with the backward trajectory clustering results 
in this paper. According to the results of quantifying backward trajectory, the precipitation events 
related to CW and ITSC accounted for a large proportion. Due to low relative humidity, the 
precipitation events generated by CW in the study area were mostly drizzle or light rain. The 
water vapor carried by CW has evaporated almost completely in the long-distance transport 
process. The proportion of precipitation brought by ITSC was the highest, and the precipitation 
events were mainly concentrated in the monsoon period (from June to September). The air mass 
from ITSC, carrying a lot of water vapor, passed through the south and southwest of China at a 
high speed and brought a lot of precipitation to the northwest inland areas. Claus et al. (2021), Wu 
et al. (2021), and Wu and Bedaso (2022) deemed that the backward trajectory tracking method 
and second-order variable of isotopes (1.e., d-excess) could reflect the degree of imbalance in the 
evaporation-condensation process of water vapor source areas; and they also thought that 
d-excess could serve as an important indicator for tracing the source of water vapor. Based on the 
results of d-excess wind rose diagrams (Fig. 6), we found that precipitation in the monsoon 
margin was mainly influenced by monsoon moisture and westerly moisture. Although we set 
sampling sites in the monsoon margin, more sampling sites for a longer duration and larger spatial 
coverage are still needed in the future. 


4.3 Influencing factors of precipitation isotopes 


Dansgaard (1964) indicated the correlation between precipitation isotopes and environmental 
factors. The air humidity and temperature are important meteorological conditions for affecting 
precipitation isotopes (Wang et al., 2017; Li et al., 2022; Rhoujjati et al., 2023). We researched 
the correlation between precipitation isotopes and meteorological factors in this study and found 
that 5!8O value was positively correlated with temperature and negatively correlated with relative 
humidity in the study area. Compared with temperature and relative humidity, precipitation had a 
smaller correlation with the precipitation isotopes due to the low intensity of precipitation in the 
study area. The sources of water vapor are also considered to be important factors affecting the 
variation of precipitation isotope component (Tian et al., 2007; Crawford et al., 2013; Baker et al., 
2015). We found that precipitation predominantly formed by dry continental water vapor tended 
to have higher d-excess value, while precipitation predominantly formed by moist marine water 
vapor tended to have lower d-excess value in the monsoon margin. In addition, the transport 
duration and distance of air mass are also vital to modify the isotopic component in precipitation. 
In this study, the relative humidity values of the air mass from CW, CNW, OA, and CP were 
62.37%, 63.58%, 70.79%, and 79.60%, respectively (Table 3). The d-excess values of the four 
sources gradually decreased with the increase of air mass relative humidity. The d-excess value of 
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the precipitation from CSE was the highest (84.90%), which corresponds to its relatively high 
temperature and short transport path with less rainout. The air mass from ITSC had the highest 
temperature when compared with air mass from continental area, and the amount of precipitation 
from ITSC was relatively considerable. The heavy isotopes in water vapor are continuously 
eliminated throughout the transport path, as a result of air mass movement. Relative humidity was 
relatively low in the air mass from PR, CNW, and CW. With the extension of transport path, the 
d-excess value of the three areas steadily decreased. 


5 Conclusions 


In this study, we analyzed the precipitation isotope data gathered in the monsoon margin for the 
entire year 2022. Based on this, we used the HYSPLIT-based backward trajectory model, PSCF 
analysis, and CWT analysis to investigate the water vapor sources in the monsoon margin. The 
LMWL of the monsoon margin was 5D=7.345'%O-1.16, defined by regression analysis. The lower 
slope and intercept compared with the GMWL (SD=8.175!8O+10.56) indicated the existence of 
strong regional evaporation in this area. Furthermore, the atmospheric water line slope of summer 
half year was higher than the slope of winter half year because of the high humidity in summer. 
We found that the effects of temperature (positive effect) and relative humidity (negative effect) 
to 8'8O value were significant through analyzing the relationship between the 5'8O value in 
precipitation and meteorological factors at daily level. The component of 5°H and 5!8O isotopes 
was also greatly influenced by the air mass trajectory. Eight areas accounted for the majority of 
water vapor sources in the study area: CNW (13.98%), CW (32.20%), CP (5.51%), CSW (4.24%), 
CSE (3.81%), ITSC (33.47%), OA (4.24%), and PR (2.54%). The precipitation in the study area 
was mostly caused by the water vapor carried by ITSC because of the high moisture content of 
ITSC. Long water vapor transport distance increased the loss of precipitation isotopes. 
Understanding the atmospheric process of precipitation is important for managing water resource 
in arid and semi-arid areas. The research time of this study is relatively short (only one year), so 
the future study should take multi-year precipitation data to analyze the linkage between the 
composition of precipitation isotopes and water vapor sources in arid areas. 
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